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ABSTRACT 

A torus of reduced differential rotation can form in the inner ^ lOpc core of active 
galactic nuclei incurring a density enhancement that can account for obscuration of X-rays 
in Seyferts when the initial inner core to black hole mass ratio > 50. The same density 
enhancement and reduction in differential rotation can also lead to dynamo growth of 
poloidal fields which attain a magnitude ~ lO^G when accreted onto the central engine. 
As radio jet models often employ poloidal fields as agents in extracting power for the 
jet luminosity, we suggest that jetted AGN might require this poloidal field production. 
Although the poloidal field would be originally produced in the obscuring torus, jetted 
objects are less likely to have obscuring tori: The poloidal field would only aid in powering 
jet emission after it accretes with the torus matter to the central engine. Thus, only during 
the relatively short torus accretion time scale could there be both a jet and and torus. 

Subject Headings: Galaxies: Active, Jets, Magnetic Fields, Radio Continuum: Galaxies 
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1. Introduction 

At the centers of galaxies, and particularly in active galactic nuclei (AGN), are likely 
black holes (BH) which power the central luminosity by accretion (e.g., Rees 1984). Time 
variability and the estimated accretion efficiency seem to require black hole masses ^ 
IO^Mq in many jetted and jet-free sources. Little is known about the detailed mass 
distribution in the "inner core" (IC) ^ lOpc regions of all galaxy types, but there in our own 
Galaxy the rotation curve seems to incur a change from Keplerian to flat (Genzel & Townes 
1987). The transition region is quasi-rigidly rotating. This leads to a density enhancement 
over a purely Keplerian curve, as we show below, and may explain the presence of the 
cicumnuclear ring (CNR) of the Galaxy and X-ray obscuring Seyfert tori (Yi et al. 1994, 
Duschl 1989) in accordance with unified models of AGN (e.g. Antonoucci, 1993). 

We show that such a region also favors the dynamo production of poloidal magnetic 
field (PMF) to a magnitude which, when accreted onto an AGN central engine, is likely 
~ lO^G — in agreement with that inferred by other equipartition estimates (Begelman et al. 
1984). Such PMF is often required in jet models (Blandford & Znajek 1977, Lovelace et al. 
1987, Appl & Camenzind 1993, Lynden-Bell 1995). The PMF dynamo growth time scale 
is much smaller than the torus depletion time, so whether significant PMF is produced 
depends only on the initial to BH mass ratio. Though the jet PMF originates in the torus, 
the field can only play a role in jets after accreting to the central engine. Only during the 
short time when the torus is depleting could there be both a jet and torus in one object. 
This is consistent with the fact that direct evidence for tori comes mainly from radio quiet 
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objects (Urry and Padovani, 1995), but more data are needed. 

2. Estimation of Time Scales and Adiabatic Black Hole Growth 

Rotation curves of spiral galaxies show quasi-rigid rotation in the inner ~ kpc, and 
flat rotation curves outside ~ kpc (Oort 1978, Binney & Tremaine 1987). Models which 
account for observed galaxy gas rotation curves seem to require (Binney & Tremaine 1987) 
(i) a central BH (ii) a stellar IC within a few-lOpc, (iii) a more diffuse nuclear bulge of 
several lOOpc, and (iv) an isothermal sphere of dark matter on kpc scales. Here we are 
interested in (i), the IC region sub-structure to the overall rotation curves, which is likely 
similar for all galaxy types. 

We flrst show that a typical BH grows by accretion slowly compared to orbital time 
scales, but rapidly compared to the IC relaxation time, so the hole's growth is nearly 
adiabatic: At any time, the accreting region will be in an approximately steady state if 
the viscous time scale r^js "C r^, where Tg is the BH growth time scale. An estimate for 

"-^v^T; ^(iF^y^'^^Qy^ Uxio^W [^^) ' ^ ^ 

where is the rotational velocity, Vr is the radial velocity, and Rd is radial extent of the 
accretion flow. We estimate the BH growth time by assuming that the objects radiate at 
Eddington luminosity Lsad- For a central BH of mass Mh, Lsdd = ^t^G{Mh /Mq)c/ k = 
1.3 X 10^^{Mh/Mq) erg/sec, where k is the Thomson opacity ~ 0.4cm2 g"^. As the AGN 
radiates at ~ Lsdd we have Msdd = 2.2 x 10~^e~^ MHy"r~^ , where e < 1 is the efficiency 
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factor for radiation by accretion. Then 

Tvis <^rg< Mn/MEdd = 4.5e xlO^r. (2) 

We must also compare r^, the gravitational relaxation time of the IC region, with 
To, the orbital time scale. When the BH is small, its radius of influence Rjn ~ GMh/cf'^ 
satisfies < Ric, where Rjc is the IC radius and a is the stellar velocity dispersion, so 
that Tr is independent of Mh- When Rj^ > Ric the relaxation time depends on Mh. In 
the Fokker-Planck approximation, for the case Rjn < Ric we have (Binney & Tremaine 
1987) 

Tr ~ lO^Vr ((T/lOOkm s-^)^(10^Mq pc'^/p,^), (3) 
where is the IC mass density. For ^ Ric the relaxation time is 

~10^V (MH/lO^M0)3((7/lOOkms-^)-3(7V/lO9)-\ (4) 
where N is the number of IC stars. 

Because of (2), (3) and (4), Tq -C tj^ -C r^., where th is the Hubble time ~ 10^*^ 
yr and Tq < 10^ yr(i?/c/10pc)(cr/100km s"-*^)"-*^, the adiabatic approach is appropriate. 
Young (1980) considers the adiabatic evolution of an isothermal sphere with a growing BH. 
Quinlan et al. (1995) confirm the results of Young (1980) and extend to non-isothermal 
spheres. As applied to the IC, these are reasonably consistent with a total mass dependence 
on radius given by 

Mtot{r) = Mh[1 + m{r/RicT], r < Ric, 
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where m = 



Mtotir) = {Mh + Mic)r/Ric, r > Rjc, 
Mjc/Mh is the ratio of BH to stellar core mass, and n ~ 3. 



(5) 



3. Gas Density 



The surface gas density is given by S(r) = 2H{r)p{r), where H{r) is the half-scale 
height of the gas and p(r) is the gas density. Following the standard treatment (e.g., 
Pringle 1981) we take the viscosity coefficient to he v = ^vtH, where we assume 7^1 
is the viscosity parameter and vt is the turbulent eddy speed. Combining this with the 
conservation of gas mass and angular momentum, the surface density satisfies (Pringle 
1981, Yi et al. 1994) 

dS/dr + /(r)E = ^(r), (6) 

where /(r) = {r^Q' /Q)' /{r^fl' /Q), g{r) = -[MQ/{27T^v^r)][2{Q/rQ') + 1], and Q is the 
angular velocity determined by the potential. The ' indicates d/dr. The solution of (6) is 



given by E(r) = exp 



"- /; dsfis)] (So + /; dsg{s)e^p [/^ dXf{X)]) , where E(ro) = Eq. 
To find the relationship between the surface density and the stellar mass density we note 
that for circular orbits, O = (GM/r"^)^/^. Plugging in for /(r) and g{r), using (5), S = 
f/RiCi and A = m-^/"'G we have 



So i?2^([n - 3] m5^+2 -352) 



( - 3] mS^^^ - 35^ ) K 



1 + 



^0 Jmi/»Sn 



(n + l)A" + l 
^^\A + A"+i)V2 

(7) 



where K = Mm5(i-Vn) (G'j^^^^^)i/2/(27r7-i;|,). 



The third term in (7) is negative, so the second term gives the correct order of mag- 
nitude. Fig. 1 shows the surface density at r = Rjc ~ Rt, the edge of torus, (outside of 



which the potential drops) for 4 values of n as a function of m and vice versa. For large m 
and n ~ 3 the density can be 100 times the Keplerian (i.e. when m <C 1) value. That n ~ 3 
is consistent with Young (1980) and Quinlan et al. (1995). This density enhancement can 
account for X-ray obscuration in Seyferts (Yi et al. 1994). 

4. Rigid Rotation and Poloidal Field Growth 

As shown above, the enhanced density torus results from reduced differential rotation. 
We can calculate the reduction in fl' by setting fl = {GMtot{f) / r^Y^'^ . Using (5) with 
n = 3 gives VL'{Ric) = -{3/2)G^/^RJ^^^mI/^ /{I + mfl'^. For m > 50 this gives a factor 
> 7 reduction in fi' from the Keplerian value. This region may be important for the 
dynamo generation of PMF for radio sources as we now describe. 

The mean field dynamo theory (Moffatt 1978, Parker 1979) splits the velocity and the 
magnetic field into mean (V, B) and fluctuating (v', b') components. The time evolution 
of mean magnetic field is given by (Moffatt 1978, Parker 1979) 

aB/at = Vx[VxB + Q!B-(AM + /5)VxB], (8) 

where a and /? are the helicity and diffusion dynamo coefficients and are functions of the 
turbulent velocity. A non- vanishing ex. is the result of buoyant eddies rising in an upwardly 
decreasing density gradient, while conserving their angular momenta. The magnetic vis- 
cosity. Am, satisfies Am <S 

For small Q', the "ct^ dynamo" (Moffatt 1978) is favored, because the maximum 
growth rate depends on o? as we will see. Sufficient reduction of Q' means that the radial 
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PMF produced by the a. effect is not slieared into toroidal field. Simulations (e.g., Donner 
& Brandenburg 1990) show that dipole modes are favored in such a dynamo, in contrast 
to the a — J] dynamo often applied to disks (Parker 1979). 

To estimate when the dynamo is favored, we work in cylindrical coordinates (r, 0, z) 
and write V = rl7(r, 2)e<ji, B = B(^{r^z)k^(^ + Bp, where P indicates the poloidal (r, 2) 
component, and indicates the toroidal component. Bp = V x A{r^z)^^. Assuming a 
and /3 are constant, (8) can be written (Moffatt 1978) 

dB^jdt = r(Bp-V)n - a(V2 - r'^)! + /^(V^ - r-^)^^, (9) 

and 

dA/dt = aB^ + /3(V^ - r'^)!. (10) 

An CK^ dynamo will dominate the a — Q dynamo when the second source term on the right of 
(12) dominates the first, that is when a/r ^ \ri}'\. Prom Parker (1979), we have a ~ OAvt- 
For a turbulently supported dust torus, observations require vt/V^ ~ 0.5 where = rfi 
(Krolik & Begelman 1988, Urry & Padovani 1995). Using these and (5), the requirement 
near r = i?/c ~ Rt becomes 0.20M^/^(1 + m)^/^ > (3/2)Mj/^(l + m)-^/^, or 

m>6.5. (11) 

When (11) is satisfied, we can ignore the first term on the right of (8) near r = Ric- 
We capture the essence of an dynamo, by considering the case when the 2;-variation 
dominates and assuming solutions of the form B(j),A oc rexp(7t + kzz). Plugging these 
into (9) and (10) gives 

= aAkl - PB^kl (12) 
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and 



aB^ - pAkl, 



(13) 



so that the growth rate is given by 



Re[7] 



(14) 



The growth rate is positive if kg < 3a/ p. The maximum growth rate is Re[7]maa; = 
(9/8)q;^//3, showing the dependence. The second term on the right of (14) provides a 
more conservative estimate. 

Let us see why an dynamo favors PMF. For an AGN torus, the fluctuation scale 
is determined by the size of dust containing clouds. The dust must be in clouds because 
it could not survive if the random velocities of > lOOkm/s were thermal. Thus the cloud 
size Tc, satisfies Vc < Rti where the torus radius Rt is the scale for variation of the mean 
quantities. We can estimate the dust cloud size from observations of the CNR of our 
Galaxy (Genzel & Townes 1987), which is thought to be similar to the AGN dust tori 
(Krolik & Begelman 1988). These observations (Genzel & Townes 1987) show clouds with 
0.1 < < 0.25pc. Now we estimate the PMF produced: Setting -Re[7] ~ k^a and using 
kzA ~ Bp, with a ~ OAvt and /? ~ {l/3)rcVT in (12) and (13) gives Bp ~ for the 
dynamo. The analogous equations to (12) and (13) for the a — Q dynamo, derived 
by keeping the first term on the right of (8) and dropping the term linear in a, give 
[Bp/B^]a-n ~ rc/Rr. Thus [Bp/B^]^^/[Bp/B^]o,-a ~ Rt/tc > 50 - 100, showing that 
the dynamo favors PMF in comparison to the a — O dynamo. This ratio is important 
for jet models, particularly when the resulting luminosity depends on B^ (e.g., Blandford 
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&; Znajek 1977). Thus a factor of 50-100 in the field corresponds to a factor of 2.5 x 10^ 
to 10"^ in the jet luminosity. 

5. Poloidal Field in AGN 

For PMF to be produced by a working torus dynamo, (11) must hold. In addition, 
the dynamo growth time must be shorter than the torus lifetime. That is, 

Td < rgMr/MH. (15) 

Note Td for the torus must be less than the field diffusion time: < (lOpc)^//? ~ 
100kpcV(100kms-i lOkpc) ~ lO^r. For a density pr ~ 5 x IQ-i^g/cm^ (Krolik & 
Begelman 1988) and radius 5pc with height 2.5pc, Mt ~ 10'^ M^. From (2) with e ~ 0.1, 
~ 5 X lO^yr. Thus violating (15) requires the extreme case of Mh > lO^My so that (8) 
is more stringent. 

The larger m is the greater the density enhancement and produced PMF. Equipartition 
between turbulent and magnetic energy densities gives an upper limit to the field. For pT ~ 
5 X 10"^^g/cm^ (Krolik & Begelman 1988), corresponding to m > 50, and vt ~ lOOkm/s 
the turbulent energy is {1/2)ptv'^ ~ 2.4 x 10~'^erg/cm'^. Setting this equal to Bp/Sn we 
have Bp ^ 8 x 10~^G. The field is accreted to the central engine as the torus depletes. The 
radial component of PMF is then sheared, but only that fraction of the field at a much lower 
scale height than that of the ~ 5 — lOpc torus. The z-component is unaffected by the shear. 
An estimate of the accreted PMF can then be made from fliix freezing. For an ion-electron 
accretion disk with height to radius ratio Hd/Rd ~ 1/50 and density Pdisk ~ 10~^ cm""^ at 
r ~ lO^'^ cm (e.g., Celotti et al. 1992), flux conservation implies that PMF will accrete to 
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Bp,disk < Bp^T[Pdisk{Hd/Rd)/pT(HT/RT)?/^ ~ 8xl0-2[10-8(l/50)/5xl0-i8(l/2)]2/3 ^ 
1.5 X lO^G, in agreement with standard estimates (Begelman et al. 1984). 

Only when the torus depletes by accretion can the field produced there move to the 
central engine and play its role in jets. As any jet formation time in the engine is likely 
much shorter than the torus depletion time, only during the latter time scale can a jetted 
object show both a jet and a torus. For a lO^M — lO^M© hole accreting at the Eddington 
rate with e ~ 0.1 and a torus mass of lO^M©, this accretion period lasts ^ 10^ — 10^ yr. 
Note also that a torus scale dynamo need not determine the final scale of the magnetic 
field, but would mediate the initial energy extraction from the rotating central source. 
The initial jet flow could drag the resulting field to kpc-Mpc scales as in Blandford &; Rees 
(1974). 

6. Conclusions 

Reduced fl' tori in the central < 10 pc regions of AGN can obscure X-rays and incur 
dynamo production of PMF likely required for jets, when m > 50 initially. PMF growth 
in the torus allows angular momentum transport, and the torus will accrete to the central 
engine carrying its field. PMF transport and torus depletion are associated processes and 
jetted objects would be less likely to show obscuring tori. Jet-free AGN would either have 
an intial m < 50 or have their jets beamed away from us. There are a few radio loud 
quasars (RLQ) or Seyferts with inferred BH masses > lO^M© (e.g. NGC 5548, Krolik et. 
al. 1991). The absence of direct evidence for dust tori in the latter may be consistent with 
our paradigm, but may require non-adiabatic analysis. 
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Figure Caption 

Figure 1: (a) Surface density at the torus' outer edge as function of m forn = 0.01, 1, 2.7 , 3 
going from the bottom to the top curves, (b) Surface density as a function of n for 
m = 0.01, 1, 10, 100 from the bottom to top curves. 



12 



